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^ -. ABSTRACT 

o 

o . 

, The increase in spin frequency as the burning atmospheres of Type I X-ray 

^ I bursts cool provides a strong constraint on the radius of the underlying neutron 

star. If the change in spin frequency is due to a change in the thickness of the 
atmosphere, the radius of the star must exceed 2>GM/(? for any equation of state 
■ and approximately 3.5GM/c^ for most physically reasonable equations of state. This 

constraint arises because the direction of the Coriolis force for radial motion reverses 
for R < GI(?{M -|- I /B?'). Furthermore, the marked change in the magnitude of 
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, the Coriolis force near compact stars provides a straightforward explanation for why 

^ \ the frequency of the quickly rotating bursters shifts by the same amount as the slow 

' rotators; they are slightly more massive, 1.6 versus 1.4 M©. 

o ■ 

O 
O 

,1^ • 1. Introduction 

O ' If material accretes onto the surface of a neutron star sufficiently slowly, a layer of fuel develops 

' and then suddenly ignites producing a burst of x-rays known as a Type I X-ray burst. These 

Ci ■ thermonuclear flashes each release about 10^^ ergs and repeat on a timescale from hours to days 

(Lewin, van Paradijs & Taam 1995, Bildsten 1998). Although only one source (SAX J1808. 4-3658) 
^ , exhibits periodic variation in its quiescent emission (Wijnands &: van der Klis 1998, Chakrabarty 

&; Morgan 1998), during the bursts themselves the emission is quasiperiodic {e.g. Strohmayer, 
Zhang & Swank 1997), apparently due to rotational modulation of the inhomogeneities in the 
thermonuclear burning. Strohmayer &; Markwardt (1999) find that the oscillations in the cooling 
tails of the X-ray bursts from 4U 1702-429 and 4U 1728-34 are nearly coherent with Q ~ 4000 and 
consistent with an increase in the modulation frequency as the burning layer cools. 

Joss (1978) calculated the first numerical models of thermonuclear burning on the surface of 
a 1.4 Mq neutron star with a radius of 6.6 km. He found that the initial surface layer of helium 
expands from a thickness of about three meters to one of thirty meters during the onset of the 
burst (also Bildsten 1995). Strohmayer et al. (1997) argue that this increase in radius is sufficient 
to account for the frequency shifts observed during the bursts due to the conservation of angular 
momentum. 



To use the conservation of angular momentum to explain the change in frequency of the 
emission as the burning region of the atmosphere (or hotspot) expands and then contracts, one 
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assumes that the specific angular momentum of a fluid element is conserved during the burst. The 

magnetic fields of the neutron stars are likely to be weak (otherwise there would be periodicities 
in the persistent emission); consequently, this is a viable assumption. However, the relationship 
between the specific angular momentum of a fluid element and its position is more complicated in 
general relativity than in Newtonian theory {e.g. Abramowicz, Miller & Stuchlik 1993). In the 
case of the Schwarzschild spacetime, the angular momentum of a fluid element is 

dL = dmOr^ = dmfi- — — (1) 

^ 1 - 2M/r ^ ' 

where is the angular velocity of the element as measured by an observer at infinity. This 
expression and its generalizations result in the reversal of the centrifugal forces {e.g. Abramowicz 
& Prasanna 1990, Abramowicz 1993) and gyroscopic precession [e.g. Nayak h Vishveshwara 1998) 
in curved spacetimes. 

Rather than using Equation 1 to obtain the change in the angular velocity of the hotspot, 
§ 2 calculates the change the angular velocity in a frame that rotates with the neutron star. 

The reasons for this treatment are two-fold. First, the calculation is quite straightforward and 
illustrative. Second, the effects of frame dragging are likely to be important for neutron stars, so 
rather than deriving a new definition for the radius of gyration to include frame dragging, frame 
dragging may simply be included in the metric. § 3 applies these results to Type I X-ray bursts, 
and finally § 4 summarizes the conclusions. 



2. The Coriolis Force in Curved Spacetime 

Evaluating the Coriolis force in the curved spacetime surrounding a neutron star is 
straightforward once the metric is specified. The forces that affect the motion of the neutron star 
atmosphere are most easily expressed in a frame that rotates with the star. In this frame the 
nonzero components of the metric tensor are 

500 = 1- — + J^Vsin2ef-3 +lj (2) 

503 = 530 = J^r^sin^ 0(^2^ + 1) (3) 
/ 2M\-i 

511 = - 1-— (4) 



r / 

522 = (5) 
533 = -r^sin^e. (6) 

where x'* = [t, r, 0, 0]. J = 7^2 which accounts for the effects of frame dragging to lowest order in 
(Landau h Lifshitz 1987), and the various quantities are given in geometric units with c = G = 1. 
For rotational frequencies , v < 400 Hz, this is a good approximation {e.g. Morsink Sz Stella 1999) 
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but for u ^ 600 Hz, the induced quadrupole moment of the star may introduce a significant but not 
dominant correction to the metric (Cook, Shapiro & Teukolsky 1994, Sibgatuhin &; Sunyaev 1998); 
this issue will be addressed in the more detailed sequel. 

The Coriolis force for material moving radially manifests itself through the connection term, 

rgi = if/'530,i + <7'°5oo,i) (7) 



l + 2//r3 M 
r^J ^ 1 - 2M/r~r 



(8) 



and T\i = Fqq = 0. The first term yields the standard "Newtonian" Coriolis force reduced by the 
effects of frame dragging, and the second term arises because the usually afferent gravitational 
acceleration has an azimuthal component for material moving radially in the rotating frame. 

The atmosphere expands with a four velocity v whose components in the rotating frame are 
found by setting u-u = \, v-v = l and u • v = ^ where u points along the time axis in the 
rotating frame (i.e. it is four velocity of an observer stationary in the rotating frame). Due to 
pressure gradients, the trajectory of gas does not obey the t and r components of the geodesic 
equations. In the non-rotating frame, these pressure gradients also have a (f) component. However, 
in the rotating frame, the gas is free to follow a geodesic in the (f) direction, yielding 

^^ = -.,.=5(i-?^-^)(:-HM)-',o(a3) (, 

where (3 is the magnitude of the three-velocity measured locally by an observer stationary in the 
rotating frame {i.e. u) and 7 = l/-\/l — 

If the radial motion is nonrelativistic and small relative to the initial radius (as for the 
expansion of the atmosphere), the result is easily integrated. After transforming the result to the 
nonrotating (distant observer's) frame, it is 

dlnJ^oo / 3M I \ f 2M\-i , , 

= -2 1- — 1-— • (10) 



dlnr V r J \ r ) 

This result in the absence of frame dragging (i.e. / — 0) could also be obtained by using the 
definition of Abramowicz, Miller h Stuchlik (1993) for the radius of gyration in general relativity. 
Comparison of this result with Equation 1 yields the following definition for the radius of gyration 
with frame dragging included to lowest order in fi, 

1M\-\^ ( Ir + M \ , , 

1-— exp --—5-^7 (11) 



^ Vl - 2M/r V r ) \ 4 W^r^ 

As the mass of the star and its moment of inertia vanish, the familiar result obtains. However, in 
the vicinity of a neutron star, the correction may be large. 
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3. Application to X-Ray Bursts 

Strohmayer &; Markwardt (1999) find that during a X-ray burst, the observed frequency of 
the modulation exponentially decreases toward a value which is constant from burst to burst for 
each source. They interpret this asymptotic frequency as the spin frequency of the star. At the 
beginning of the burst, they argue that the atmosphere expands and slows due to the Coriolis 
force, and as it cools and deflates the spin frequency increases again through the Coriolis force. 
This simple inverse correlation between the spin frequency and the thickness of the atmosphere 
restricts 

R>ZM + —r (12) 

for the neutron stars associated with 4U 1636-54, KS 1731-26, Aql X-1, the galactic center source, 
4U 1702-43 and 4U 1728-34 ( R and M are the mass and radius of the neutron star). For a 
uniform density sphere in the Newtonian limit I = ^MR'^. Ravenhall &; Pethick (1994) find that 
the moment of inertia for the FPS equation of state (Friedman & Pandharipande 1981) is well fit 

by 

I = 0.21Mi?2 (i - —) . (13) 



R J 

In the context of this equation of state, R > 3.49M. 

Gumming & Bildsten (2000) have examined the structure of the atmosphere of the X-ray 
burst in further detail. By understanding the magnitude of the change in the frequency of X-ray 
bursts as well as its sign, tighter constraints on the radius of the neutron star may be obtained. 
The the proper thickness of the layer (A^i.gg) as calculated by Gumming & Bildsten (2000) for 
gs = 1.89 X 10^^ cms~^ and the coordinate thickness (Ar) are related by 

Ar 7 1. R f 2M\ , , 

- = 2.1xlO-WA.,.,-(l-— ) (14) 

where the variation in the surface gravity with mass and radius is included and 
2.1 X 10~^cm~^ = ggC^- The change in observed frequency Au for u = 300 Hz is 

A. = -0.252 Hz^^^ fl - — - ^) (15) 
20 m 300 Hz M V R R^ J ^ ^ 

It is also straightforward to add the relativistic corrections to the more detailed results of 
Gumming & Bildsten (2000) (C&B) which yields 



AOoo 



R r R 



GR M \ M 



C&B 



-1 



3M I \ An 



^ R R^J n 



(16) 

C&B 



where the effects of changing surface gravity are also included by assuming that the scale 
height of the atmosphere is inversely proportional to the gravitational acceleration; Gumming &; 
Bildsten (2000) assume that R/M fa 4.82. 
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Figure 1 shows the corrections to the calculated values of AJl including both the effects 
of general relativity and varying the surface gravity. If the effects of the curved spacetime are 
included, the result is given by the middle curve, and if frame dragging is also included the lower 
curve traces the result. The expected change in the observed angular velocity as the atmosphere 
expands vanishes for R = 3M (without frame dragging) and R = 3.49M (with frame dragging). 

Table 1 of Gumming & Bildsten (2000) presents the properties of the oscillations from 
Type I bursts. Specifically, the sources fall into two types those with u ^ 300 Hz and those with 
u 600 Hz. Both sets of objects typically have Au ~ 1 Hz. Miller (1999) argues that those 
objects with faster oscillations simply have two burning fronts active simultaneously, so the actual 
rotational frequency of the star is ~ 300 Hz and the kinematics of the expanding atmosphere is 
similar to those with faster oscillations. An alternative possibly is that the 600-Hz oscillations 
correspond to a 600-Hz rotational frequency and that Au/v is smaller due to a different value of 
M/R. 

To wit, assuming that the 300-Hz objects are 1.4 M© neutron stars with the FPS equation of 
state (Friedman & Pandharipande 1981, Ravenhall &; Pethick 1994) they have a radius of 10.75 km 
and M/R = 0.196. Using this data to normalize Equation 16 for Av/u = 6 x 10~^ gives 
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IT 



= 13.5 X 10"^ (17) 

C&B 



which is about twice as large as the value that Gumming &; Bildsten (2000) find for their 
Newtonian models. 

Using this value. Equation 16 can be used to infer a value of M/R for the fast rotators of 
0.234, modestly different from that of the slow rotators (0.246 if frame dragging is neglected). 
Without the general relativistic corrections, M/R for the fast rotators would be 0.391 which is 
more compact than is achieved for any standard equation of state. Using the value of M/R which 
includes frame dragging yields M = I.6M0 for the fast rotators and R = 10.3 km. 

If more detailed models argue that the normalization (Equation 17) is unphysically large, the 
conclusion would be that the objects have R ^ 3M. Specifically, if the modeling of the evolving 
atmosphere by Gumming & Bildsten (2000) is precisely correct, R = 8.1M for the slowly rotating 
bursters and R = 5.7M for the fast rotators. If the mass of the neutron stars exceeds a solar mass, 
this requires a very hard equation of state such as the TI and MF results (Pandharipande, Pines 
& Smith 1976). 



4. Conclusions 

The Coriolis force in the curved spacetime near a neutron star is substantially smaller than 
the Newtonian value and reverses its direction as the radius of the star approaches 3GM/c^. For 
typical neutron-star parameters of i? = 10 km and M = 1.4 Mq, the fully relativistic Coriolis force 
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Fig. 1. — The rclativistic corrections to the expected spin-down from Type I X-ray bursts. The 
upper slanted line gives the Newtonian result, and the cross shows that it is normalized for a 
neutron star with R = 10 km and M = 1.4 Mq. The middle line gives the rclativistic correction 
without frame dragging, and the lower curve traces the correction including frame dragging. 
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is only 30% of the Newtonian value. This contrasts with the assertion of Strohmayer (1999) that 

the correction to the Coriolis force is the same order as the gravitational redshift, i.e. 20 - 30% 
due to time dilation alone. Although time dilation does contribute, it does not dominate. Time 
dilation alone would not predict, the change in the sign of the Coriolis force for r ~ 3M. 

The bulk of the effect arises from the fact that the specific angular momentum of a fluid 
element has a minimum at the photon circular orbit, which has r = 3M for a non-rotating 
star. This also manifests itself through the reversal of the centrifugal force also at r = 3M 
{e.g. Abramowicz Sz Prasanna 1990, Abramowicz 1990, Abramowicz 1993, Sonego &; Massar 1996) 
and the fact that compact stars may have moments of inertia which exceed |Mi2^, the Newtonian 
value for a sphere with constant density {e.g. Hartle 1967, Chandrasekhar Sz Miller 1974). 

Because the Coriolis force sensitively depends on the ratio of the mass to the radius of a 

relativistic star. Type I X-ray bursts provide a strong constraint on the equation of state of the 
underlying neutron star. Specifically, since the oscillation frequency of the burst increases as the 
atmosphere cools and presumably shrinks, the radius of the star must exceed 3 GM/c^. If frame 
dragging is included, the constraint is stronger, R > 3.49 GM/c^. Furthermore, this sensitivity 
provides a simple explanation for the fact that the quickly rotating bursters, 4U 1636-54, 
KS 1731-26, Aql X-1 and the galactic center source, spin up by the same amount as the slowly 
rotating bursters, 4U 1702-43 and 4U 1728-34. The fast rotators are simply slightly more massive, 
1.6 M0 versus 1.4 Mq. Further studies of the thermonuclear burning in the atmospheres of Type I 
X-ray bursts will make them a precise probe of the spacetime geometry surrounding rotating 
neutron stars. 

General relativity presents gravity as an inertial force; therefore, it is not surprising that 
Newtonian notions of inertial forces do not apply in strong gravitational fields. Gravity strongly 
affects the Coriolis force which is important in the evolution of Type I X-Ray bursts; consequently, 
these bursts provide a strong constraint on the spacetime geometry surrounding accreting neutron 
stars. 
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like to thank Lars Bildstcn, Greg Ushomirsky, Alessandra Buonanno, Lior Burko, Scott Hughes, 
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